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THE  BRAGG  CRYSTAL  SPECTROMETER  AS  A DIAGNOSTIC  TOOL 


1.  INTRODUCTION 

The  Bragg  crystal  spectrometer,  consisting  principally  of  a diffraction  crystal  and  a 
detector  sensitive  to  the  diffracted  radiation,  is  useful  for  determining  source  characteris- 
tics in  the  X-ray  region  of  the  spectrum.  Not  only  can  it  be  used  for  identifying  spectral 
features,  but  also  for  determining  quantitatively  the  flux  in  spectral  lines  and  in  the  con- 
tinuum. Because  its  response  is  sensitive  to  the  source  polarization  direction,  it  can  be 
used,  with  some  limitations,  as  a polarimeter.  If  sufficient  information  on  the  diffraction 
characteristics  is  known,  it  can  also  be  used  to  determine  spectral  line  shapes. 

In  the  following,  we  will  .first  outline  the  theory  of  crystal  diffraction  phenomena. 

We  will  then  develop  a mathematical  description  of  the  interaction  of  source,  crystal,  and 
detector.  Using  this  description,  we  will  discuss  employment  of  the  Bragg  crystal  spec- 
trometer for  determination  of  the  properties  of  source  spectra  and,  simultaneously,  define 
the  required  spectrometer  characteristics.  Finally,  methods  for  determining  the  spectrome- 
ter characteristics  will  be  given  with  examples  of  results  obtained. 


2.  CHARACTERISTICS  OF  CRYSTAL  DIFFRACTION,  THEORY 
2.1  Bragg’s  Law 

Consider  a small  parallelepiped  of  a crystal  lattice.  This  parallelepiped  is  defined  as 
a unit  cell  containing  lattice  points  only  at  its  comers.  If  one  of  the  corners  is  chosen  as 
the  origin,  the  location  of  an  adjacent  comer  (lattice  point)  in  the  unit  cell  is  given  by 
one  of  the  vectors  A,  B,  or  C.  The  volume  of  the  unit  cell  is  then  given  by 

Y=(AXB)-C.  (1) 

The  location  of  all  other  lattice  points  in  a perfect  crystalline  block  of  such  unit  cells 
(referenced  to  a corner)  may  be  specified  by  the  vector, 

r = nA  A + nBB  + r^C,  (2) 

where  the  n are  integers  of  the  form 

"A  =0.1 Na  -1 

nB  = •••>  Nfi  - 1 

Sp  ” 0, 1, ...,  Njn  — 1. 

The  total  number  of  lattice  points  in  the  block  is  then  N = 
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To  describe  the  diffraction  (scattering)  phenomena,  we  introduce  the  vector, 

S = s - s0,  (3) 

where  sQ  and  s are  unit  vectors  in  the  direction  of  the  incident  and  diffracted  beams, 
respectively.  For  such  a three-dimensional  array  of  scattering  centers  (assuming  negligible 
absorption  and  rescattering  and  unit  index  of  refraction),  the  intensity  of  the  diffracted 
beam  will  be  essentially  that  which  we  would  expect  from  a three-dimensional  grating. 
That  is  intensity  I will  be  proportional  to  [1] 


where  X is  the  wavelength  of  the  diffracted  beam. 

For  Na,  Nb,  and  Nc  » 1,  the  intensity  of  the  diffracted  beam  will  be  strongly 
peaked  and  is  given  by 

I«(NaNbNc)2  = N2.  (5) 

The  conditions  necessary  for  these  maxima  to  occur  are 

S-A  =h\ 

S‘B  = k\  (6) 

SC  = ex, 

where  h,  k,  and  e are  integers  (known  as  Miller  indices).  Now  S is  a vector  normal  to  the 
reflecting  plane  and  of  magnitude  2 sin  0 (0  being  one-half  the  scattering  angle);  thus,  the 
conditions  of  Eq.  (6)  become 


2A  sin  0 cos  a = h\ 


2 B sin  0 cos  p = ftX  (7) 

2 C sin  0 cos  7 - CX, 


where  a,  P,  and  7 are  the  angles  between  the  normal  to  the  reflecting  plane  and  the  vec- 
tors A,  B,  and  C,  respectively.  The  direction  of  S is  therefore  uniquely  defined  for  given 
h,  k,  and  V,  and,  further,  there  is  a constant  “spacing”  of  the  lattice  planes  given  by 


A B „ C 

DhkU  = -r  cos  o = — cos  P = — cos  7 


(8) 
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so  that  our  maxima  occur  for 


X = 2 Dhk%  sin  0,  (9) 

which  is  known  as  Bragg’s  law  of  diffraction.  If  h , k,  and  C have  a common  factor,  we 
arbitrarily  remove  it,  so  that 


dh/m,k/m,t/m  ~ mDhkZ> 


(10) 


and  Bragg’s  law  becomes 


mX=2dhlm,klm,Zlm  sin9n 


(ID 


which  we  may  interpret  as  mth-order  diffraction  from  the  plane  (him,  k/m,  C/m)  with  the 
appropriate  2d-spacing. 

For  the  particular  case  of  cubic  crystals  of  one  element, 

Dhki  = A(h2  +fe2  + C2)'1'2  (12) 

since  A,  B,  and  C are  orthogonal  and  of  the  same  magnitude.  For  a more  general 
parallelepiped  unit  cell,  Dhk g and  the  reflecting  plane  normal  are  most  easily  determined 
by  using  the  reciprocal  lattice.  We  define  vectors 


a = (B  X C)/V 
b = (C  X A)/V 
c = (A  X B)/V 


so  that 


a*  A = b*B  = c'C  = 1. 


Then  the  vector 


r„  = ha  + kb  + Cc 

is  in  the  direction  of  the  normal  to  the  reflecting  plane  ( h , k,  C)  and,  further. 


Q3) 


(14) 


(15) 


|r„  | - (16) 

In  practice,  for  parallelepiped  crystals,  it  is  therefore  possible  to  specify  an  approximate 
spacing  d and  to  obtain  a slab  of  crystalline  material  cut  and  ground  so  that  the  surface 
of  the  slab  is  parallel  to  the  reflecting  planes. 
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Because  the  index  of  refraction  to  xrays  for  most  materials  is  less  than  unity,  i.e., 
n = 1 - 6 (where  5 is  small),  Bragg’s  law,  as  expressed  in  Eqs.  (9)  and  (11),  is  not  com- 
pletely accurate.  For  wavelengths  considerably  less  than  any  critical  absorption  wave- 
lengths [2] , 


2nmec2mH 

* 1.36  X 10~6pX2,  (17) 

where  wavelength  X is  in  angstroms  and  crystal  density  p is  in  g/cm3.  The  deviation  of 
the  index  of  refraction  from  unity  is  usually  small,  but  for  very  accurate  wavelength 
determinations,  we  must  modify  Bragg’s  law  so  that  our  dispersion  relation  becomes 

m\  = 2d(l  - 5/sin2  9m  ) sin  0m  (18) 

where  6m  is  the  angle  measured  in  vacuum.  The  error  made  by  using  Bragg’s  law  can  be 
expressed  as 


— = 5.44  X 10~6  pd2,  (19) 

* sin2  0m 

which  for  most  materials  is  ^ 10~3,  far  from  absorption  edges  (regions  of  anomolous 
dispersion).  Nevertheless,  if  the  wavelengths  are  to  be  determined  very  accurately,  the 
modified  Eq.  (18)  must  be  used  in  lieu  of  Bragg’s  law  as  expressed  in  Eqs.  (9)  and  (11). 


2.2  The  Small  Cry  stall 

The  above  equations  may  be  derived  solely  by  the  use  of  physical  optics.  They  are 
by  no  means  complete,  because  we  also  will  require  a measure  of  the  diffracted  intensity. 
We  introduce  the  differential  scattering  cross  section  do/df2,  defined  by 


do  _ Energy  scattered/time/solid  angle 
dSl  Energy  incident/time/area 


The  scattering  of  xrays  is  due  almost  completely  to  electrons,  so  we  consider  scattering 
by  an  assemblage  of  electrons  located  instantaneously  at  positions  x;.  For  the  case  of 
our  small  crystalline  block  (negligible  absorption  and  rescattering  and  unity  index  of 
refraction),  the  scattering  is  then  given  by  [3] 


(n 

u 


da 

cm 


= <|£  — e(2»r'/X)s-K;|2)sin2  @ 
j me<? 


(20) 
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where  0 is  the  angle  between  the  observation  direction  and  the  electric  field  vector  of 
the  incident  radiation.  The  average  is  taken  over  all  possible  values  of  x;  (a  time  average). 
(This  average  will  not  be  important  to  our  treatment  here,  but  such  averaging  is  needed  in 
analysis  of  temperature  effects.)  For  the  case  of  the  electric  field  perpendicular  to  the 
diffraction  plane  Ea , sin2  0 = 1;  for  the  case  of  the  electric  field  lying  in  the  diffraction 
plane  E” , sin2  0 = cos2  2 0m,  dm  being  the  Bragg  angle  of  diffraction.  For  now  we  con- 
sider only  the  former  case  E° . 


We  now  reference  the  electron  positions  to  the  center  of  the  lattice  points  of  Sec. 
2.1;  then  the  position  of  each  electron  is,  from  Eq.  (2), 

x = nA  A + nBB  + ncC  + R,  (21) 

where  R is  the  electron  position  relative  to  the  lattice  point  or  scattering  center  (atom)  to 
which  it  belongs.  Then  the  summation  in  Eq.  (20)  becomes 


£ 


e<l  c(2m7\)S-x.  - g2 
o 2 

mec * mec 


^e(2Jri/X)nAS-A^’ 
nA  nB 


e(2m7X)nflS-B 


£ 

nc 


,(2m/X)nCs-c  ^ g(2m7X)S-Rfc 


(22) 


where  we  have  assumed  that  the  crystalline  block  consists  of  atoms  of  one  kind  only  and 
where  Rfe  is  the  position  vector  of  the  feth  electron  in  that  atom  (lattice  point)  designated 
by  nA(=  0,  1,  ....  Na  - 1),  nB(=  0,  1,  ...,  NB  - 1),  and  nc(=  0,  1,  ...,  Nc  - 1).  Since  the 
summations  are  independent,  we  may  treat  each  separately.  Now, 


„ J2tu7X)Nas-a  , 

y g(2m7X)nAs-A  =£ ± 

nA=0  e(277i/X)S-A  _ 1 

e(ni/\)NA S-A  sin  ^ ^ S.A  j 
g(ffj7X)s-A  sin  ^1.  S.A  j 

and  similarly  for  the  summations  over  nR  and  nc.  The  summation  over  k may  be 
replaced  by  an  integral  of  the  form 


(23) 


f = 


Jp(R)e(27r,/X)s'R  d3/?, 


(24) 


where  p(R)  is  the  number  density  of  electrons  in  an  atom  relative  to  its  center.  The 
quantity  f is  called  the  atomic  form  factor  or  atomic  scattering  factor. 
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By  inserting  Eqs.  (23)  and  (24)  into  Eq.  (20),  we  find,  for  a block  consisting  of 
atoms  of  one  kind  only,  that 


ft 


\mec 


2 |fl2  f S-A),m2(w.  f S-b)  sin2(wc  * S-c) 

™2  (l  S'A)  s,n2(f  S‘B)  Sin2(ls'c) 


(25) 


If  more  than' one  kind  of  an  atom  is  contained  in  our  block,  we  use,  in  lieu  of  f in  Eq. 
(25),  the  structure  factor,  which  follows  immediately  from  the  above  treatment  and  Sec. 
2.1, 


/ye(2Tr,/X)s*ry , (26) 

j 

where  the  sum  is  over  all  the  atoms  in  the  unit  cell.  Vector  r;  is  defined  as 

Tj  = UA]A  + UBj  B + ucj  c,  (27) 

where  UAj;  UBj,  and  UCj  (<  1)  are  the  distances  along  A,  B,  and  C,  respectively,  in  the 
unit  cell  where  atom  j is  found.  For  diffraction  by  the  plane  ( h , fc,  K),  Eq.  (26)  becomes 

Fhkt  = ^f/*i(hVASkUB^UV\  (28) 

i 

Now,  the  total  power  in  the  diffracted  beam  is  just 

p°  = J indn. 

Since  the  diffracted  beam  is  extremely  narrow,  we  let 

dn  = dyd& 


and 


S = + 5S, 


where  5S  is  small  and 


S/j'A  = h\ 
SH  • B * feX 
S„-C  * CX. 


6 


This  result  is  for  a particular  incident  angle  9.  By  using  the  reciprocal  lattice  vectors  a,  b, 
and  c from  Eq.  (13)  and  letting 


6S  = X(£aa+£6b+gcc),  (30) 

we  find, 

^5S-A  = t rga,  (31) 

and  similarly  for  the  other  dot  products.  We  now  wish  to  integrate  the  power  over  all 
incident  angles.  That  is,  we  wish  to  find 

Jpff(0)d0. 

Now  the  element  dydfi  in  Eq.  (29)  is  effectively  an  angular  area  element  in  a plane 
perpendicular  to  s.  We  may  express  small  changes  in  s produced  by  small  changes  in  y 
and  j3  as 

5s  = byey  + Spe0, 

where  is  a unit  vector  perpendicular  to  the  diffraction  plane  and  e lies  in  the  diffrac- 
tion plane.  (See  Fig.  1.)  We  note  that  s-8s  = 0.  Now,  from  Fig.  1 we  see  that  a small 
displacement  in  0 produces  a change  in  Sg  of 

5s0  = - 5 9 cos  2 9ey  + 6 9 sin  20ea, 

where  ea  = s. 

Thus,  since  S = s - s0, 

5S  = (5y  - 59  cos  29  )ey  + 5(3ep  + 59  sin  2 0ey. 

So,  the  volume  element  in  our  Cartesian  coordinate  system  ( ey , fy,  ea)  is  given  by 

dydflda  = dyd$d9  sin  29.  (32) 

Now  5S  is  also  given  by  Eq.  (30),  so  a differential  volume  element  is  also  given  by 

[(AadgJ  X (Xbd#ft)]  •(\cdgc)  = X3(a  X b)-c dgadghdgc.  (33) 
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Fig.  1 — The  relationship  between  the  (e  e^,  eQ)  coordi- 
nate system  and  the  incident  and  diffracted  direction 
vectors.  For  appreciable  intensity  in  the  diffracted  beam, 
0 =5  0Q.  Note  that  S = sin  (0  + + cos  + 90^®a’ 


Since  diffraction  takes  place  only  over  a small  range  of  angles,  we  let  F(20)  become 
F(29q)  and  sin  20  become  sin  20o,  where  Oq  is  the  diffraction  maximum.  From  Eqs. 
(31),  (32),  and  (33),  the  integral  of  the  power  over  all  incident  angles  is  given  by 


J*  P°(6)dd  =fA 


|F(20O)|2X3 


(a  X b)‘C 

sin  20  o 


X 


sin2(AT4  irgA ) 
sin2(7rga) 


sin2  (NBngb)  sin2  (Ncirgc) 
sin2(7 Tgb)  sin2(7rgc) 


(34) 


For  Na,  Nb,  and  Nc  » 1,  the  sin2  terms  in  the  denominator  become  their  argument  and 
the  integrals  may  be  extended  to  ±°°.  Equation  (34)  then  becomes 


I 


P°(0)d0  = 1°. 


!F(20o)|2X3 


naNbnc 
V sin20q 


(35) 


If  n is  the  number  of  unit  cells  per  unit  volume,  we  find  that 


|F(20o)|2X3 


n2(vol) 

sin20o 


(36) 


where  (vol)  is  the  volume  of  our  crystalline  block.  This  discussion  is  found  to  apply  very 
well  to  “powder  diffraction,”  where  many  small  crystals  are  employed. 


2.3  The  Mosaic  Crystal 

We  now  consider  a large  crystal  composed  of  small  blocks  of  the  type  discussed  in 
Sec.  2.2.  The  blocks  are  assumed  optically  decoupled  from  one  another,  and  the  entire 
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crystal  is  assumed  to  be  so  thick  that  no  radiation  is  transmitted  through  (mosaic  crystal). 
In  this  case,  absorption  of  incident  and  diffracted  xrays  must  be  taken  into  account.  We 
therefore  introduce  an  attenuation  coefficient  into  the  intensity  relation: 


/Sr - is  (-f-)  mo0 


|2*3 


n2A 


2p  sin20o 


(37) 


where  p is  the  linear  absorption  coefficient  and  A is  the  cross-sectional  area  of  the 
incident  beam,  and  we  have  used  Eq.  (36).  Now  the  total  energy  per  time  incident  to 
our  crystal  is  Then,  for  the  special  case  of  the  mosaic  crystal, 


X3n2 
2p  sin20o 


(38) 


Note  that  if  we  accumulate  diffracted  radiation  E,  as  a crystal  is  rotated  at  constant 
angular  velocity  co,  we  find  that 


E = j P(0)dt  = w"1  Jp(0)dd. 


(39) 


The  ratio  Eu/Iq  is  known  as  the  integrated  reflection  coefficient,  integrated  reflectivity, 
or,  more  loosely,  reflectivity  from  a crystal  face  R0. 


Inserting  now  the  polarization  dependence,  from  Eq.  (20)  and  the  following  discus- 
sion, we  find  that  for  the  mosaic  crystal, 


and 


*0 


\F(20o)\2 


X3n2 
2/j  sin  20  o 


(40a) 


R*  = E u 
0 ^ 

Jo 


Rq  cos22 Oq. 


(40b) 


These  equations  are  incomplete,  however,  since  they  consider  neither  secondary  nor 
primary  extinction.  Both  forms  of  extinction  occur  for  angles  of  incidence  near  the 
Bragg  angle;  i.e.,  they  are  diffraction  processes.  Primary  extinction  takes  place  within  a 
small  block  and  is  due  to  the  reduction  of  incident  beam  intensity  at  a given  layer  in  the 
block  by  previous  diffraction  of  other  layers  and  by  interference  between  the  incident 
beam  and  the  rescattered  diffracted  beam.  Secondary  extinction  takes  place  because  the 
beam  incident  to  a small  block  has  been  depleted  by  diffraction  of  other  (previous)  blocks; 
interference  does  not  occur  because  the  small  blocks  are  optically  decoupled  in  a mosaic 
crystal. 
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2.4  The  Perfect  Crystal 

Diffraction  from  a large  perfect  crystal  must  properly  take  into  account  absorption, 
index  of  refraction,  and  rescattering  of  the  diffracted  beam.  To  do  this,  we  depart  from 
the  above  procedure*  and  consider  a plane  sheet  of  similar  atoms  regularly  spaced.  We 
let  a plane-parallel  incident  beam  strike  the  plane  at  an  angle  0.  The  beam  is  then  dif- 
fracted from  the  plane  at  essentially  the  same  angle. 

We  let  the  plane  of  atoms  be  the  x - y plane;  the  position  of  each  atom  in  the  plane 
is  then  given  by  a position  vector,  so  that 


r = xev  + yev . (41) 

We  let  the  y - z plane  be  the  diffraction  plane.  Then  the  incident  beam  is  described  by  a 
wave  vector,  such  that 

k = ~~  (cos Oe  -sinOe^).  (42) 

a y 


Incident  to  each  atom,  then,  we  have  wave  radiation  of  the  form 

A0ei^‘-K',)=  AQei^t-Kycos8'> 


(43) 


where  to  = 2irc/\  and  /10  is  the  amplitude  of  the  incident  beam.  Each  seatterer  (atom) 
will  scatter  the  incident  beam  in  all  directions  in  the  form  of  a spherical  wave.  The 
amplitude  of  the  scattered  wave  will  be  (a-radiation  only,  i.e.,  E°  is  perpendicular  to  the 
diffraction  plane) 


- Aof(20 ) 


^i(uit-KycosB-Kp) 

P 


(44) 


at  a distance  p from  the  seatterer.  If  we  now  observe  the  scattered  radiation  at  a point 
in  the  diffraction  plane  at  an  angle  0 above  the  plane  of  scatterers  and  a distance  R from 
the  origin  (i.e.,  the  position  vector  of  the  observation  point  is  R = R (cos 0ey  + sinfle, )) 
we  find  that 


p = [R2  + r2  - 2Ry  cosO] 1/2. 


(45) 


Now  R » r and  for  our  purposes  the  factor  p in  the  denominator  of  Eq.  (44)  becomes 
simply  R.  This  simple  substitution  may  not  be  made  in  the  argument  of  the  exponential 
term,  however.  For  this  purpose  we  expand  Eq.  (45),  finding  (to  second  order  in  r and 

y) 


p R - y cos 0 + r2(l  - sin20  cos2d)/2E 


(46) 


♦The  procedure  we  employ  is  similar  to  that  of  R.  W.  James,  The  Optical  Principles  of  the  Diffraction  of 
X-Rays,  Vol.  II  of  The  Crystalline  State  Sir  L.  Bragg,  ed.,  pp.  52-66,  Cornell  University  Press,  Ithaca, 
New  York,  1965. 
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where  we  have  substituted  r sin <p  for  y.  The  amplitude  of  the  scattered  beam  from  one 
scatterer,  given  by  Eq.  (44),  is  now 

„2  ei[wf-Kfi-Kr2(l-sin20cos20)/2i?] 

-Aof(20)  (47) 

mecl  n 


We  define 

amplitude  of  the  scattered  beam  from  a single  scatterer 

Q — . — ■—  ■ — — - ■—  # ( 4 C 

s amplitude  of  the  incident  beam 

We  then  see  that  the  total  contribution  at  our  observation  point  due  to  all  the  scatterers 
in  the  x - y plane  is 


27T  r00 


-»/ / 

*'f\  */A 


q rdrd<p, 


(49) 


where  N is  the  area  density  of  scatterers  and  we  have  assumed  the  plane  sheet  of  scatter- 
ers to  be  large.  From  Eqs.  (47)  and  (49),  we  find  that 


q=-Nf(20) 


e2  ei(w/-Ki?)  r2n  r°° 


mec~  " -'0  ■’0 

The  integrals  are  quite  straightforward  and  yield 


f f £~(iKr2/2R)(l-sin20cos20)  (5Q) 

•'A  •'A 


q = iNf(20)  e 

KM  /I  “ 


i(co(-kR) 


k sind 

which  is  again  a plane  wave.  This  equation  can  be  rewritten  as 


q = indf(20 ) 

mec2 


_ ei(wf-«R) 
sinO 


(51) 


where  we  have  replaced  k by  2ir/\  and  N by  nd  (the  product  of  the  number  density  of 
scatterers  and  the  distance  between  planes  of  scatterers).  If  we  translate  back  to  the 
origin  (R  = 0)  and  ignore  the  time  factor  (an  unnecessary  complication  in  the  following 
treatment),  we  have  finally,  for  the  ratio  of  diffracted  to  incident  amplitudes  (including 
the  effect  of  a phase  lag  of  very  nearly  n/2,  which  takes  place  on  reflection). 


- iq  = i(nd\/s\n0)f(20) 


(52) 
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Fig.  2 — Relationship  between  the  reflected,  incident, 
and  transmitted  amplitudes  and  the  reflection  planes 


The  same  sort  of  development  of  the  ratio  of  the  amplitudes  of  the  forward  scattered  and 
the  incident  beams  gives 


- iq0  = i(McfX/sinf? )/"( 0) 


(53) 


Let  us  now  consider  a crystal  composed  of  many  layers  of  such  planes,  strictly  parallel, 
and  with  all  scatterers  in  register.  Let  Tk  be  the  incident  amplitude  to  the  loth  plane  and 
Sk  be  the  total  amplitude  of  the  scattered  waves  from  all  planes  t > k.  The  fraction 
transmitted  by  each  plane  is  (1  - iq 0);  the  fraction  reflected  is  - iq.  Then,  as  shown  in 
Fig.  2, 

Sk=-  iqTk  + (1  - iq0)e^i2n,^sindSk+l  (54) 


and 


Tk  = (1  - iq0)e-(i2vlX)d&ind  Tk  - Iqe-<,27r/X>2dsin0S*+1  (55) 


where  we  have  included  phase  factors  to  account  for  the  distance  d between  planes 
(assuming  plane-parallel  radiation),  and  q , rescattering  from  the  reverse  side  of  plane  k. 
These  two  equations  lead  to  an  equation  for  the  transmission  amplitudes  only  and 
another  for  the  scattered  amplitudes  only.  Then  setting 

Tfe+l  = xTk  ancl  sk+i  =ysk 

(where  |jc|  and  lyl  are  assumed  to  be  very  nearly,  hut  less  than,  unity)  we  find, 

-e-(i2nl\)dsin0  + _ jqQ )2e-(i2»r/X)dsin0  + ?(/2rr/X Wsinfl  (5 g) 

and  an  identical  equation  for  y.  From  our  assumptions,  x = y,  and  therefore  also  S,.+l  = 
xSk.  If  we  assume  that  0 is  nearly  the  Bragg  angle  for  the  wavelength  of  the  incident 
radiation,  the  phase  factor  becomes 

c-(i27r/A)dsin0  _ p-imn-iv  ^57) 


(1  ~ iQq) 


(H- 
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where  v is  small  so  that 


e~iv  « 1 -to. 


Also,  since  \x\  is  nearly  unity, 

x * (1  - £)e-'m7r. 

Substituting  for  x in  Eq.  (56),  we  find 

£2  « qq  - (<J0  + v )2 


(58) 


(59) 


(60) 


Using  Sfe+1  = xSk  and  Tfc+1  = xTk  yields 

S0  = - iqT0  + (1  - iq0)e-i(m7r+v)xS0. 

Thus  the  ratio  of  diffracted  amplitude  S0  and  incident  amplitude  T0  is,  from  Eqs.  (58), 
(59),  and  (60), 

£o  = ~ “?  — ^ ~q  ■ - . (61) 

T0  1 - (1  - iq0 )(1  - ii/)[l  ±y/qq-  (<f0  + ^ <?0  + v ± y/(q0  + v)2  - Qq 

'rhe  ratio  of  diffracted  to  incident  intensity  is  then  just  1S0/T0I2,  and  we  choose  the  sign 
of  the  radical  to  guarantee  that  the  reflection  coefficient  is 

R°(v)  = |S0/T0|^  <1.  (62) 

We  now  consider  the  special  case  of  negligible  absorption.  If  the  crystal  does  not 
absorb,  q = q and  q0  and  q are  real.  Then, 

= ~ <7  (63) 

T0  q0  + v ±v/(90  + u)2  " qi 


and  for  |q0  + i>l  < q, 

R°(v)  = 1.  <64) 

Again,  outside  of  this  region  the  sign  of  the  radical  is  chosen  to  ensure  the  validity  of 
Eq.  (62).  Note  that  the  peak  of  Ra(v)  occurs  for  v = - q0.  If  0o  is  the  Bragg  angle,  the 
peak  occurs  for  an  angle  0o  + A0,  where 


AO 


nd\\f(0)\(e2  IrrifC2) 

- cos00  sin()0 
A 


sinf>0  cos Oq 
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corresponding  to  the  modification  of  Bragg’s  law  expressed  in  Eq.  (18).  The  width  is 
approximately 


A v = - 2 q. 


This  corresponds  to  an  angle 


50  = 2A0  -1-  . 


90 


Then,  for  o-polarization. 


50°  = 2A0 


A 20) 


AO) 


(66) 


and  for  7r-polarization,  the  angular  width  is 


50*  = 2 AO 


A 20) 


AO) 


|cos20 ) 


(67) 


For  negligible  absorption,  the  integrated  reflection  coefficients  are  then  given  by 

j*R(0)dO  = | 50, 


where  the  factor  4/3  comes  from  the  additional  contribution  of  Kq.  (63)  in  the  wings. 
Then  we  see  that 


(0)d0  = 


8 nX2 


3n  sin  20  0 


! A 20  o ) I 


(68a) 


(0)d0  = — 


8 n.y 


377  sin20o 


IA20O) 


|cos20 


O' 


(68b) 


these  are  Darwin’s  formulas  for  the  integrated  reflection  coefficients  of  a perfect  crystal 

[41. 


The  Prins  modification  of  Darwin’s  treatment  includes  absorption  by  allowing  5,  q0. 
q,  and  q to  be  complex  [5].  Then 


6 -*■  5 + j'(3 
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where 

|3  = X/H/47T 

(p  being  the  linear  absorption  coefficient).  Now 


and  we  can  set 


- (6  + m = 


<?  0 


sin  Or 


% d 

X 


(69) 


and 


- (A  + iB) 


q sinO 

T* 


(70) 


Setting  AO  = 0 - 0O  (small),  we  see  that 


cosOq  sin0oA0 


v sin0o 

f7 


Then,  by  substituting  Eqs.  (69),  (70),  and  (71)  into  Eq.  (61),  we  find 


A + iB 


T0  C ±s/c^Ia^B)(A^B) 


(71) 


(72) 


where  C = cos 0o  sin0oA0  - 6 - t'(3.  For  the  case  of  a composite  crystal  containing  more 
than  one^kind  of  atom,  f is  replaced  by  F (Eq.  (28))  in  q0,  q,  and  q , so  that  5,  /3,  .4,  B, 
A,  and  B are  also  modified. 


For  a nonpolar  crystal  q = q and  A + iB  = A + iB.  Thus 


R°o  = fRO 


( 0)d0  = - 


25 


sin20r 


F(20o) 

F(20) 


J 


1 + iB/A 


drj 


(73a) 


(where  the  polarization  factor  in  A and  B has  been  set  equal  to  unity),  and 

Rq  = R'£  Icos20o| 


(73b) 
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where  Rtf 
Icos20o|. 


is  of  the  same  form  as  R°  but  A and  B now  include  the  polarization  factor 


Some  actual  crystals  seem  to  be  described  quite  well  by  the  perfect  crystal  descrio 
tion  given  here,  especially  when  Eqs.  (62),  (72),  and  (73)  are  app/LJ  * P' 


2.5  Discussion  of  Theory 

It  is  not  our  intent  to  vigorously  derive  the  theories  of  crystal  diffraction  nor  to 

r^  d7e7haSPeCt  °f  Crystal  diffraction;  for  these  and  other  discussions  the  reader  is 
referred  to  the  appropriate  literature,*  from  which  much  of  the  above  material  has  been 

Constructing  an  x-ray  spectrometer  instrument  that  takes  advantage  of  the  disoersio, 

hi  Hv°n  laWVS  qUite  SimPle‘  We  are  horded.  Potentially,  a Seat  dej  of  flex! 

y m e choice  of  crystal  and  of  atomic  plane  spacing,  so  the  wavelength  ranee  of 
interest  can  also  be  matched,  within  limits.  It  is  therefore  possible 
lines  are  present  and  the  ratio  of  line  radiation  to  continuum  in  the  neighborhood  of  the 


dpfin,n77rtUn7ely;T  h7C  at  1<?ast  two  markedly  different  theories  (mosaic  and  perfect) 
defining  the  integrated  reflection  coefficient,  only  one  of  which  (perfect)  defines  a 

tfonToTdiffrlcti  R !hthUS  lmfOSS‘ble  t0  determine  the  sPectrUm  of  the  incident  radia- 
and  knn7«dtffr  l + th®?iy  UI)leSS  °ne  kn°WS  which  <if  either>  idealized  theory  pertains 
^d  knows  i .a  structure  factors  (including  temperature  effects,  not  included  here)  One 
must  therefore  determine  for  each  crystal  those  parameters  required. 

In  the  following,  we  explore  the  expected  results  obtained  by  a Bragg  crystal  spec- 
trometer for  various  types  of  sources  and  ways  of  experimentally  determining  the  parame- 
ters needed  to  define  the  incident  source  spectrum.  running  me  parame- 


3.  INTERACTION  OF  SOURCE,  CRYSTAL,  AND  DETECTOR 

The  typical  arrangement  of  a Bragg  crystal  spectrometer  is  shown  in  Fig  3 The 

Ta^lTo0/  !!et:ryst7are;rtd  7 an  angie  0 to  the  incident  beam;  «»* 

is  at  an  angle  2 6 to  the  incident  beam  and  views  the  diffracted  radiation  from  the  crystal. 

tion  ^Vhe°Z  dfCr'bed  earlier’  dlffractlon  from  a c^stal  depends  on  the  polariza- 
Hnn  77  rad,atlon-  We  will  see  that  in  some  respects  this  is  a complies- 

he  ; eneVT  SS  * P°tentially  useful  for  determining  the  state  of  polarization  of 

Uth  f f 10  3ny  CESe’  WC  mUSt  fmd  S°me  means  for  treatinK  it-  A further  complication 
is  that  few  x-ray  sources  are  monochromatic,  and  it  is  nearly  impossible  to  obtain  one 


D.  Van  N,.  * « 
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Fig.  3 — Typical  arrangement  of  crystal  and 
detector  of  a Bragg  crystal  spectrometer 


that  is  not  divergent.  In  the  following  treatment,  we  will  find  that  the  characteristics  of 
our  spectrometer  that  must  be  determined  in  the  laboratory  are  (a)  the  quantum  efficiency 
of  the  detector;  (b)  the  integrated  reflection  coefficient;  (c)  the  reflectivity  curve; 

(d)  the  polarization  dependence;  (e)  the  atomic  plane  spacing;  and  (f),  for  very  accurate 
identification  of  spectral  features,  the  index  of  refraction.  The  means  by  which  the 
quantum  efficiency  of  a detector  may  be  determined  is  similar  to  that  dealt  with  at  some 
length  in  Meekins  et  al.  [6]  and  will  not  be  further  considered  here;  we  will  assume  the 
detector  efficiency  to  be  known.  In  addition,  atomic  plane  spacing  d and  index  of 
refraction  may  be  obtained  to  sufficient  accuracy  by  least-squares  fitting  techniques 
applied  to  known  laboratory  spectra.  Thus,  they  also  will  not  be  further  considered;  we 
will  assume  them  to  be  known. 

In  the  following,  we  will  demonstrate  the  need  for  the  above  spectrometer  charac- 
teristics (or  lack  of  them)  when  dealing  with  specific  sources  as  well  as  the  methods  for 
their  determination  in  the  laboratory.  We  first  devise  a means  for  treating  rotations  when 
we  have  a polarized  source. 


3.1  Polarization  and  Rotations 

Consider  a source  emitting  linearly  polarized  radiation.  We  may  describe  the  radia- 
tion using  amplitudes,  which  we  take  as  parallel  (II)  and  perpendicular  11)  polarization  in 
some  reference  frame.  We  thus  describe  the  sources  as 


If  we  simply  view  the  source  with  a detector  having  unit  efficiency  to  both  II  and  1 
components,  we  in  fact  find  an  intensity,  flux,  or  power,  which  is 
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p =&*&=  Cf  + Cf . (75) 

Suppose  we  have  a detector  sensitive  to  linear  polarization,  oriented  so  that  the 
source  and  detector  reference  frames  are  coaligned.  A workable  representation  of  the 
detector  response  is  then 


(76) 


where  the  e are  detector  efficiencies.  Then  the  intensity  as  seen  by  the  detector  is 

Pm**fia~C\el+CleL  (77) 

Also,  in  our  formalism,  we  will  require  the  description  of  rotations.  Since  amplitudes 
transform  as  vectors,  a rotation  of  \p  about  the  direction  of  propagation  is  described  as, 


/ cos\p  0 sirup' 
7 sirup  cosi p , 


(78) 


where  we  must  determine  0 and  7.  We  have  set  the  coefficients  of  the  cos\p  terms  equal 
to  1 so  that-ty^  = 1 for  \p  = 0.  Now  v/^,  must  be  unitary,  so  that 

I/3I2  = I7I2  = 1 

and 

7*  = -0. 


Now,  suppose  we  perform  a rotation  between  the  detector  and  source.  Then, 

P'*t*l**j9r 


or 

P = e||fcj|  cos2 1 p + C2  sin2 0 + C1C2(0*  + (3)  sirup  cosi//]  + eJC2  sin2 + C2  cos 2\p 
- Cj  C2(0*  + 0)  sirup  cosi//] . 

From  the  Law  of  Malus,  we  know  that  the  dependence  of  intensity  on  angle  of  rotation 
1 p is  of  the  form  sin2i p and  cos2i p,  not  sirup  cosi p.  So  /3*  = -0,  and  we  choose  0 = i so 
that  counterclockwise  rotations  may  be  considered  positive.  We  also  find  that  circular 
polarization  may  be  treated  by  allowing  \p  to  be  a function  of  time.  X-ray  measurements 
are  usually  made  over  times  that  are  long  in  comparison  to  the  cycle  time  of  1 p,  however, 
and  therefore  time  averages  of  cos2\p  and  sin2i//  must  be  taken,  leaving  us  with  (Cj^  = 0), 

P=\  (f||+ei)Cf 

i.e„  no  information  on  the  polarization. 
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We  now  insert  a polarizer  (e.g.,  crystal),  which  in  its  own  reference  frame  is  described 
by 


1. e.,  diagonal,  since  it  cannot  mix  polarization  states  if  the  source  polarization  reference 
frame  is  the  same  as  that  of  the  polarizer.  If  we  now  use  a source,  a polarizer,  and  a 
detector  (which  is  rotated  by  i//  with  respect  to  the  polarization  reference  frame  of  the 
polarizer  and  source),  we  find  that  the  intensity  is 

P=.3t  tx  ty/t  / r/^.xj, 

and  we  conclude  that  A||  and  A±  have  the  same  phase,  which  we  set  equal  to  zero  so  that 
A ||  and  Aj_  are  real. 

Normally,  the  x-ray  detector  itself  is  not  sensitive  to  polarization,  so  that 

fi  = e^.  (80) 

Then,  for  the  general  case  of  a multicomponent  system  with  rotations  i//lt  2»  ^ „ 

between  polarizers  and  source,  we  would  have 

P = e[Xn^U^n  (81) 

In  the  following,  we  will  consider  two  configurations:  1.  A source,  a crystal 
(rotated),  and  a detector  (Fig.  4);  for  this  case, 

P(e^)  =e  Mtf/^12 

= eAfl(C2  cos2^  + Cf  sin2^) 

+ eA2(C^  sin2 + C2  cos2 <^)  (82) 

2.  A source,  two  crystals  (the  second  rotated),  and  a detector  (Fig.  5);  in  this  case, 

P(l,ei'1/)  = e\.XB«4l.xAm2 

= eB2(A2C2  cos 2\p  + A2C2  sin2^) 

+ eB2(A2C^  sin2 1// + A2C2  cos2 1/*)  (83) 

where  we  have  indicated  the  rotations  of  each  crystal  by  and  parameters  of  the  second 
crystal  by  Bj  and  B^. 

We  see  that  the  treatment  of  polarization  effects  and  the  interaction  of  rotations 
about  the  incident  beam  are  straightforward  with  the  use  of  these  matrix  operators. 
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DETECTOR 

K2®  / 


Fig.  4 — Diagram  of  the  polarization  reference  frame 
and  a rotated  crystal.  The  crystal  is  represented  by 
the  reflection  plane.  The  diffraction  plane  contains 
both  the  incident  and  diffracted  x-ray  beams  and  is 
normal  to  the  reflection  plane. 


3.2  The  Source,  Crystal,  and  Detector:  Operator  Approach 

We  will  describe  the  source  as  seen  by  our  detector  or  crystal  by  a function  SXtS  in 
units  of  photons  per  time  per  wavelength  interval  per  (angle)2.  The  principal  direction 
will  be  along  the  x-axis,  so  divergence  angle  vector  6 lies  in  the  y - z plane.  We  insist 
that  the  centroid  of  the  incident  radiation  and  the  x-axis  be  coincident,  so  that 


I 


SSXt5dSydSz  = 0. 


(84) 


As  noted  above,  diffraction  from  a crystal  depends  on  the  state  of  polarization  of  the 
incident  beam  and  its  orientation.  Whenever  this  must  be  considered  we  will  write 


s\tS  ~s\tS  +S{/S- 

Source  function  S is  related  to  the  amplitudes  of  Sec.  3.1  by 


(85) 


(86) 
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In  nearly  all  cases,  it  is  possible  to  restrict  the  divergence  of  the  source  to  very  small 
angles,  so  that  we  may  describe  a specific  ray  incident  on  our  crystal  or  detector  by  a unit 
vector: 


n 


ev  + 


Syev 


+ 6zez. 


(87) 


Over  a small  range  of  incident  angles,  the  quantum  efficiency  of  most  x-ray  detectors  is 
constant  and  independent  of  polarization.  Thus,  if  we  view  an  x-ray  source  with  a detec- 
tor, we  would  observe 


= J c(X)SXlSdMbyd8z 


(88) 


in  units  of  photons  per  time.  Thus,  we  may  describe  the  detector  with  a detector  opera- 
tor identical  for  both  polarization  states: 


E 


J 


c(\)d\dq^dqz 


(89) 


where  dq%dqz  indicates  integration  over  all  incident  angles. 

From  the  earlier  discussions,  we  know  that  diffraction  from  a crystal  depends  on  the 
angle  between  the  atomic  planes  and  the  incident  x-ray  beam,  and  that  the  peak  intensity 
of  the  diffracted  beam  occurs  when  Bragg’s  law  (Eq.  (11))  (or  small  departures  therefrom) 
holds  true,  i.e.,  when 


m\  ^ 2d  sinOm . 


Even  with  an  incident  beam  with  no  divergence,  there  is  a small  spray  of  radiation 
diffracted.  There  is  also  a small  range  of  incident  angles  over  which  diffraction  takes 
place  (albeit  not  at  peak  intensity).  Referring  to  Fig.  6,  we  let  the  principal  direction 
(centroid)  be  along  ex;  the  normal  to  the  atomic  planes  is  then  represented  by 


N = cosaex  + sinaev  (90) 

so  that  ez  lies  in  the  atomic  (reflection)  plane  and  normal  to  the  principal  directions  of 
both  incident  and  diffracted  beams.  If  we  represent  the  incident  beams  by  Eq.  (87),  the 
total  diffracted  intensity  over  all  diffracted  angles  will  be  proportional  to  a function  (to 
first  order  in  8 8Z,  and  0 - q>)Q(0C\)  - s e - 8y),  where  0<X)  = sin (mX/2d)  and 

= 7T/2  - a.  Referring  to  Fig.  7,  we  let  the  principal  direction  of  the  diffracted  beam  be 
er.  Due  to  mirroring  at  the  crystal,  the  small  divergence  angles  5V,  hz  become  - 8V,  - 8Z 
upon  diffraction.  Then  to  describe  the  ray  of  the  diffracted  beam  we  use 


n„  = (r?v  - 8y)ey  + (qz  - 8z)ez  + 


(t?v  - 8y)2  + (rjz  - 8Z)~ 


(91) 
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Fig.  6 — Incident  and  diffracted  beams  of  a 
crystal,  showing  the  crystal  normal  and  the  dif- 
fracted spray  about  the  centroid. 


Fig.  7 — Incident  and  diffracted  reference  frames. 
The  divergence  angle  becomes  negative  on  reflec- 
tion. 


where  rjy  and  tjz  are  the  spray  angles  in  the  corresponding  directions  for  a nondivergent 
incident  beam.  Then,  to  first  order  in  angle,  the  spread  about  the  centroid  may  be 
described  by  the  function  QI(r?v  - 5v)ev  + (tj*  - 5^)e;].  We  arbitrarily  employ  the 
normalization  and  centroid  conditions, 


(92-al 
(92-h) 

I 


Q(x)dx  = 1 


y 

f.rQt.t 


x)dx  - 0 
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J Qiyey  + zez  )dydz  = 1 (92-c ) 

and 

/(ye,  + ze2  )Q(ye,  + zez  )dydz  = 0.  (92-d ) 

Then  our  entire  diffraction  function  is 

Q[0(X) 6y] -,R0(X)-Q[(t?,  - 8y )ey  +(tjz  -8z)ez], 

where  fl0(X)  is  a normalization  factor  and  is  a local  constant. 

Suppose  we  have  a spectrometer  having  diffraction  characteristics  described  by  our 
diffraction  function  and  a detector  of  unit  efficiency.  We  irradiate  the  crystal  with  a 
monochromatic,  nondivergent  source  of  the  form 

s\tS  =S,8(X  -\0)5(6z)8(6y). 

Then  the  observed  count  rate  will  be 


n 


t 


\0)6(6z)5(6y)d6yd6z 


’Q(0(\)  - $ - 6,  ).R0(X)Q[(r?,  - 5V  )ev  + (tjz  - 6Z  )e;  | 

•dXdfr?,  - Sy)d(7)z  - 5Z) 

= S(R0(\q  )Q(0(Xq  ) - ifi). 

If  we  now  rotate  the  crystal  at  a uniform  rate  d^/dt  = oj,  we  find  that  the  total  number 
of  photons  contained  in  the  diffracted  line  is 

n — co  S,*0(X0). 

Then 


(93) 

which  we  may  use  as  an  experimental  definition  of  the  integrated  reflection  coefficient, 
in  agreement  with  Eqs.  (40),  (68),  and  (73).  From  our  above  discussions,  we  know  that 
the  diffraction  will  depend  on  the  state  of  polarization  of  the  source  and  its  orientation 
with  respect  to  the  diffraction  plane.  We  therefore  define  a crystal  reflectivity  operator 

3S 

<KCT(X,*>)  =J  dqzdqyQ»lOCK)-<f,-qy]R°(\) 

X Qa[(7v  “V®.v  + -«,>•*]  (94) 
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where  we  have  used  dqzdqy  to  indicate  integration  over  the  angles  of  the  incident  radia- 
tion. There  is  a similar  expression  for  the  crystal  reflectivity  operator  to  be  used  with 
7r-polarization.  We  have  then  a general  crystal  reflectivity  operator, 

^(X,  = d?a(X,  ^f>)  + bln(\,  ifi)  (95) 

each  part  acting  on  only  that  part  of  the  incident  radiation  that  has  the  indicated 
polarization.  We  have,  then,  established  a connection  between  our  reflectivity  operators 
and  the  polarizer  parameters  squared,  A\  and  A\.  (See  Eq.  (79)  and  the  following). 

Suppose  we  examine  a monochromatic,  distant  point  source  of  radiation,  linearly 
polarized  at  orientation  k, 

s\t8  =s;  5(X  - )8(5y ) (96) 

with  a Bragg  crystal  spectrometer.  From  such  a source,  we  will  observe  a count  rate  in 
our  detector  of 

nt  = £<R(X,  'P)S^ls . 

If  we  let  the  direction  between  the  k -polarization  incident  and  the  normal  to  the  dif- 
fraction plane  be  i//  (see  Fig.  4),  from  Eq.  (82)  we  find  that 

nt  = fa 6(\  - \0)8(8z)8(8y)d5zd8y 

•{cos2^Qc’tO(X)-vj-5ylftg(X)Q«|(7?v  -6y)ey  + (Vz  - 
+ sin2^Q7r(0(\)-^-5y]«g(X)Q^[(T?v  -6v)ev  + (^-«*)ej} 

•e(\)d\d(qy  - 8y  )d(rjz  -8Z). 

Carrying  out  the  indicated  operations  yields 

n,  =S'‘e(X0){cos2^Q(,[tf(X0)  -^1«£(X0)  + sin21//Q7T((?(X0)  - v?]fl£(X0 )} . 

If  we  now  orient  our  spectrometer  so  that  \p  = 0, 

< =S?e(Xo)Qa(0(Xo)-^]Kg(Xo) 

and  if  ^ = w/2, 

nf  = 5?e(X olQ^lflfAo)  (100) 

If  we  scan  the  spectrum  at  an  angular  rate  d<fi/dt  = w,  the  total  number  of  counts 
(photons  observed  by  the  detector)  observed  under  the  line  would  be,  from  Eq.  (99), 

n°  =S^e(X0)R"(X0)w-1  ,101) 
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m 

and,  from  Eq.  (100), 

nn  =S?e(A0)«£(A0)w-1. 

(102) 

If  we  observed  the  source  directly  with  our  detector,  the  count  rate  would  be 

nt  = S?e(X0). 

(103) 

Thus,  the  integrated  reflection  coefficient  (at  X0)  to  o-radiation  is,  from  Eq.  (101), 

= n°“>lnt 

(104) 

and  that  to  ^-radiation,  from  Eq.  (102),  is 

*S<X0)  = nnco/nt. 

(105) 

For^a  randomly  polarized  source,  we  let  Sf  - S°  + Sf,  where  S°  = Sf  in  magnitude. 

nt  = e(X0)(S?  + Sf) 

(106) 

and  the  total  number  of  counts  from  a randomly  polarized  source  is 

| nr  = €(X0)co-1[SfRg(X0)+S}rRg(X0)]. 

(107) 

From  Eqs.  (106)  and  (107),  the  integrated  reflection  coefficient  for  incident  random 
polarization  is 

Rr0(X0)  = wn^/n,  = ~ l«o<xo>  + Ko(X0)1  • 

(108) 

For  an  arbitrary  state  of  polarization  of  the  incident  radiation,  we  let 

Sf  = sj  +s|  (orthogonal). 

(109) 

Then,  if  we  perform  all  the  indicated  operations  (letting  i//  be  the  angle  between  1-polari- 
zation and  the  normal  to  the  diffracted  plane  (see  Fig.  4 and  Eq.  (82))  we  find  the 
count  rate  to  be 

nt  =f^o){QoI0(^o)-^]«g(Xo)[S<icos2|/,  +sf  sin2^  1 

+ Q7T(0(X0)-^]R5(X0)(s|lcos2t//  + Sj  sin2  \!/] .} 

(110) 

Scanning  through  the  line,  we  find  that  the  total  observed  counts  would  be 

n = [«o<XoHS/  ™s2^  + S\  sin2^)  + RJ(X0)(S«  cos2^  + S,1  sin2*)l 

(111) 
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while 

nt  = e(X0)(S}+s|'). 

Thus  we  see  that  we  may  set 

n,  = ( n° 1 + nf11)  cos2i p + (n^11  + nj^)  sin2i£ 
n = (nai  + nnlt)  cos2i //  + (rtal1  + rtni-)  sin2i^ 
and 

~nt  = nj  + rif  (112) 

which  will  simplify  our  mathematics  as  we  continue  to  more  complex  sources. 

The  application  of  our  operators  to  idealized  line  sources  is  reasonably  straightfor- 
ward even  when  polarization  effects  are  included.  We  have  found  experimental  defini- 
tions for  the  integrated  reflection  coefficients  for  idealized  sources  (Eqs.  (104),  (105), 
and  (108)),  which  may  be  compared  with  the  theoretical  definitions  in  Eqs.  (40),  (68), 
and  (73). 


4.  APPLICATIONS  TO  REAL  SOURCES  AND  USE 
AS  A DIAGNOSTIC  TOOL 

In  this  section,  we  will  apply  the  Bragg  crystal  spectrometer  to  realistic  sources 
using  our  operators.  Simultaneously,  we  will  find  those  spectrometer  characteristics 
required  for  determining  as  much  as  possible  about  the  source  spectrum.  We  have 
already  mentioned  the  use  of  Bragg’s  law  (Eq.  (11)),  from  which  spectral  features  as 
observed  by  our  spectrometer  may  be  identified.  Quantitative  aspects  of  the  spectrum 
are  somewhat  more  difficult  to  determine,  however,  as  we  shall  see. 

Real  sources  are  divergent  and  nonmonochromatic  and  may  also  be  arbitrarily 
polarized.  We  will  describe  these  sources  by 

r*“sx,*+sJ/a-  <113> 

Then  count  rate  nn( 1 is 

n?1  = js^6d6id5vQo[0(X)-^-6y]R«(X)Q<,[(r?y  - 6V  )ev  + (r?* 

• e(X)dXd(i7y  - 6y)d(rj;  - bz).  (114) 
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Suppose  the  divergence  is  not  too  large.  Then,  since  R0  and  e are  both  local  constants 
(reasonably  far  from  absorption  edges  in  both  crystal  and  detector)  and  the  crystal 
diffracts  in  only  a small  range  of  angles,  Eq.  (114)  becomes 


n°tl  = Rg(Xo)e(Xo)Js^5d62d5ydXQ[0(X)  ].  (115) 

If  we  let  the  source  emit  both  line  and  continuum  radiation,  where  the  continuum  is 
locally  constant  in  wavelength,  then 

= sjjg (line)  + (continuum),  (116) 

and  it  follows  that 


n°L  = -R^Xo)e(M{/ s{^d62d8yd\Q°l0a)  - * - hy  ) 


f sX?6dM^JV[0(*)-*>-6y]dX 


(117) 


Now  from  Bragg’s  law  for  first-order  diffraction, 

dX  = 2d  cos OdO 


and,  since  is  very  nearly  equal  to  0 , 


*£(WX0>  Jsjfod8,d6ydXQ[0(X)-*-6y]  +2dcos*S^J 


(118) 


I f1 

where  Sjj  is  the  number  of  photons  per  wavelength  interval  per  time  contained  in  the 
continuum.  If  we  scan  over  the  line  at  a rate  co  = d»p/dt,  (the  elapsed  time  being  At)  we 
find,  just  as  for  Eq.  (101), 

noi  =K“(X0)e(X0)(co-lS,u  + At  2dcos^,)  (119) 

where  S jL  is  the  total  number  of  photons  per  time  contained  in  the  line  over  which  the 
scan  took  place.  While  the  identification  of  a particular  line  in  a spectrum  (including  its 
diffraction  in  several  orders)  is  usually  unambiguous,  identification  of  the  continuum  is 
often  not,  so  we  include  all  the  diffraction  orders  by  setting 

naL  = co-1R«(X0)e(X0)S‘/'  + At  2d  cos  * ^ ^ R"(X0/m)c(X0/m)S^/m  , (120) 

m 
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where  m is  the  order  of  diffraction.  Far  from  a line,  but  not  so  far  that  R0,  e,  <p,  and 
SjJ’  are  appreciably  different, 


n°v  = 2d  cos * £ ^ R°0(\0/m)e(\0/m)S{colm  (121) 

m 

Then  the  number  of  counts  in  the  line  is 

nolL  = + continuum)  - Atn^fcontinuum).  (122) 

Usually,  a linear  interpolation  of  nt  (continuum)  is  performed  to  determine  the  con- 
tinuum contribution  to  the  spectrum  under  the  diffracted  line.  If  so  determined,  At  nt 
is  much  more  nearly  the  true  contribution  under  the  line.  Now  the  total  counts  in  the_ 
line  observed  by  the  detector  will  be,  from  Eqs.  (112),  (120),  121),  and  (122), 

nL(line)  = w-le(X0H[Kg(Xo)S}L  + «S(X0  1 coS2^ 

+ [Rg(X0)sjli  + «S(X0)S^]  sin2  *}  (123) 

which  is  identical  to  the  result  obtained  from  a monochromatic  distant  point  source. 

(See  Eq.  (111).)  The  continuum  count  rate  in  the  spectrum  is 

nf  (continuum)  = 2d  cos *>  £ ± e(X0/m){[fl£(X0/m)Sj£/m  , + /{J(X0/m)Sg/m  ,)  cos2^ 

m 

+ [«0<X0 t + Rno^ol"')Sfolm  ,1  sin**  }.  (124) 


If  we  have  reason  to  believe  that  the  source  is  randomly  polarized, 

SU  = SjL  m 1 SL  and  SLC  =SIIC  = -I  Sfx, 

so  that, 

nL  =co-lf(X0)S/,'K^(X0)  (125) 

and 

nf  = 2d  cosy?  ^ f(Xo/",)s>x0/mKo(Xo/m'  (126) 

m 

where  we  have  used  Eq.  (108).  If  no  higher  order  continuum  is  present  (or,  if  any,  such 
has  been  removed). 

nf  = 2d  cosipe(X0)^o/?Q(X0).  (127) 
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Even  for  a source  with  unknown  orientation  \p  with  respect  to  the  plane  of  diffraction, 
it  is  possible  to  find  the  total  contribution  to  the  spectrum  from  the  sum  of  both 
polarizations  just  as  if  the  source  is  randomly  polarized.  We  make  two  spectral  scans, 
one  at  some  orientation  4/  and  the  other  at  i //  ± n 12.  Then  the  total  of  the  two  scans  of 
the  continuum  becomes 


£nf  = n?ic  +n?llc  + n*ic  +nf"c 
or 

H n?  = 4d  cos^e(X0)/f"0(X0)S^0,  (128) 


where  we  have  again  removed  any  higher  order  contributions,  and  Sf  , = sic,  + S-!|c 
Similarly,  for  the  line  radiation,  0 A°'  A°r 


= 2fl'(X0)SfoJ-le(X0).  (129) 


Thus,  to  find  the  source  spectrum  due  to  the  total  of  both  polarizations  from  the  Bragg 
crystal  spectrometer  results,  we  must  know  the  integrated  reflection  coefficient  Rr^\)  due 
to  a randomly  polarized  source.  If  we  had  instead  taken  the  difference  between  the  two 
scans, 

Anf  = 2d  cos*  cos  2^  £ ± e(X0/m)[fi^(X0/m)  - ^(X0/m)]  (S^/m  , - S^/m  , , 


(130) 

wAn'-  = e(X0)[*g(X0) -K*(X0)][S,iL -s“L]  cos  2^  (131) 


where 


A nf  = nf 1 


^ i n/2 


and 


AnL 


nL ^ 12 


At  this  point,  we  note  that  angle  for  the  continuum  need  not  be  the  same  as  the 
corresponding  angle  for  the  line(s).  We  may  now  write 


A ^ = I^(^o)-«o^q)1 
ZnL  2 K'(X0) 


(S 


ll. 

t 


s 


l. 

I 


cos  2 , 


(132) 
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and  if  we  were  certain  that  no  higher  order  continuum  were  present  (or  that  any  such 
continuum  had  been  removed),  y 


if , «;<»„»  _ <,  , 


Inf 


2Rr0a0) 


oC 

b\0, 


cos2 1 p . 


(133) 


By  producing  a number  of  such  ratios  for  various  angles  $ + £ (£  varied)  we  can  deter 

m,„e  * (An  - 0 for  f - l„ 4 - t odd).  W,  furtor  „ou>  tot  if  „ - 

Brewster  s angle),  the  theories  on  crystal  diffraction  (Sec.  2)  give  R™  = 0.  Then  (X0  = 2d 


Ant  <S&-S£> 

~c  = ^ 9 cos2^ 

^0 1 


(134) 


for  *he.  ra‘ios  of  the  integrated  line  intensity.  For  this  particular  instance 
u A u ’ the  P°larizatlon  ratio  AS/2S  may  be  determined  without  any  knowledge 
about  the  crystal.  The  more  general  case  (arbitrary  *>),  however,  requires  Rr,  (X0)  and 
0'A0'  ~ Ko(Ao)>  > e.,  any  two  of  the  integrated  reflection  coefficients  (from  Eq.  (108)). 

We  include  one  other  case  here.  Some  sources  are  so  small  in  angular  extent  that 
they  may  be  considered  to  be  of  the  form 


SX,6  =Sx,S(5,)6(6y). 


(135) 


We  will  consider  here  only  line  radiation,  since  the  contribution  due  to  continuum  may 
be  readily  removed.  Then  y 


SX,6  = (^f  + sj[f)5(6j5(5y).  ( 

Let  us  now  assume  that  the  wavelength  distribution  in  the  line  of  our  source  mav  be 
given  by,  * 


and 


= T(X)SA 
Sjjf  = T(X)sfL 


where  T(X)  is  a positive  definite  function  normalized  to  unity.  Then  from  Eos 
and  (115),  ’ M ' 

nt  «/*(*)!  cos2«A  + {S^t^^+S^c/*^)]  sin2 

where  we  have  let 


(137) 
(112) 

(138) 


u°(v)  =Jr(X)fi"(X)Q°[0(X)  - ^]dx 

^7r(^)=/7’(X)f?J(X)Qff[0(X)-^ldX. 


(139) 
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i ' V'~t  ^ \Y) 

where  V = ((jo  + ipr)/2  T,  ,.  (14° 

I(X)  is  quite  br°ad,  its  determinatiorf  d^pends^  S°™etirnes  very  useful,  but  unless 
then  we  must  solve  an  integral  equation  ' k"0'vledge  of  Q’  and  Q«.  Even 

“ "ne’ " Wh“h  — »»)  will  be  Ute'C  ofTith:^L"’ay  ^ 
extremely  usefu^d^^'m'S'usTfu ll11*, .th?  Bra8e  'TstaJ  spectrometer  is  an 

'Z2ST' * \ ■>  I-  by  omVck  of 

available,  given  limited  knowledge  of  th*  cn . * te  a hlferarchy  of  source  information 

saar.1 

(m  one  case  with  very  little  calihratin  source  polarization  mav  bp 

integrated  reflection  coefficients  to  nr»i  • J1  L^e  sPectrometer)  provided  that  th*> 

2^5-  <?«>  * derepnt^^SiSrr  ‘“"V  ~ boib,e 
separate  nearby  Imes,  given  the  reflectivity  ££££?•■  * ‘° 

DETERMINATION  of  the  crystal  parameters 

5.1  Absolute  Integral  Reflection  Coefficients 

Suppose  we  describe  the  souree  be  used  ,or  calibratio„  by 
SM  ' <Sjf  * S,‘l )S (\-\o). S1C  , glc 

""  the  derector  viewing  diffracred  ration  from  , erysra,  Jg,  obse„e 


(141) 


nt  = £<R(\,  ^>)s 
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or,  from  Eqs.  (82),  (89),  and  (94), 

n,  = ]'{[S^5(X-Xo)+S^6]Qa[0(X)-^-6y]Kg(\) 

• Q°  [toy  - 6y  )ev  + (x]z  - sz  )e2  ] + [s"^5(X  - X0)  + s{f6  ] 

• ^[^(X)  - 6y  JJR5(X)e[(i7y  - 5y  )ey  + to2  -62)e2j} 

• e(X)dXd(rj2  - 62)d(r?y  - 6y  )d62d6y  (142) 

where  we  have  set  i//  = 0. 

Assuming  that  the  are  constants  locally  independent  of  wavelength  and  scanning 
the  spectrometer  in  $ at  rate  co,  we  find  the  total  count  in  the  line  and  continuum  to  be 

(e*'*  = 1), 

n(l)L+c  = e(XQ)  Ra0{\0)(u-1  j S^d62d8v  + 2dAf  cos*  j Sj&d6,d5y) 

+ Rg(X0)(w_1  j sjg" d62d6y  +2 dAt  cos * Jsj[f6dS2d5yj  (143) 

where  At  is  the  time  required  to  make  the  spectral  scan.  If  we  now  rotate  the  spectrome- 
ter by  7r/2  about  the  incident  beam,  we  similarly  find  (e1^  = i); 

n(i)L+c  = e(X0)  *g(X0)  (cu-1  J S^d8zd&y  + 2dAt  cos*  | sj[fgd62d6yj 

+ fig(X0)  (co'1  ( sjgdS2d8y  + 2d  At  cos * j Sj^d52d5v)  . (144) 


It  will  become  apparent  later  that  including  both  the  continuum  and  the  line  creates 
a needless  complication.  We  might  use  the  continuum  only,  but  the  number  of  photons 
diffracted  in  some  of  the  following  procedures  will  be  very  small.  The  continuum  can  be 
determined  by  interpolation  between  measurements  made  on  each  side  of  the  diffracted 
line.  The  diffracted  continuum  produces  count  rates  of 

n,(l)c  = e(X0)2dcos*  Kg(X0)  | S^6d52d5y  + R£(X0)  |sj[f6d52d6y  (145) 

and 

nt(i)c  = e(X0)2dcos*  Rg(X0)  | sj[f6d52d5y  + Rj(X0)  | sjfsd6tdSy  . (146) 
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We  now  let 


which  becomes 


h = 


n(i)MC  - n(i)c 
n(l)L+C  -ntlf 


+ Ro  fsisdS2dSy 
KfSisd8^d8y+BoJS!sd8^d8y 


(147) 


(148) 


Now  we  examine  the  spectrum  with  a second  crystal  (parameters  £01 , Qj , and  Qj ).  Then 
the  observed  count  rates,  from  Eq.  (142),  are  (i p = 0) 


nf(l)L+c  =e(X0)^1  | Js^Q®[0(Xo)-v? -Sy]d5zdSy  + 2d  cos*  J s{fsd6zd8, 
+ Rn01  | JsJfQ*[0(Xo)  - * - ?>y]dbzd5y  + 2d  cos*  f S^sd8zd6 


(149) 


and  on  either  side  of  the  line, 


n((l)c  =e(X0)2dco^f/?2i  | SxfsdM5.v  + Boi  | S\rid8*d8y)  • (150) 


Keeping  this  same  crystal  in  the  first  position,  we  examine  the  diffracted  beam  with  our 
original  crystal  (Fig.  5).  Then, 


nt  = £6((X,  x)^i(^. 

or,  from  Eqs.  (89)  and  (94)  (ignoring  polarization  effects), 

nt  = |sA/5d6id6><?1[0(X)-^-6vl£oi(x)®il(7v-6v)ev  + 0?*-S,)qJ 

• d(ny  -Sy)d(t}2  -6,)Q[0(X)-x-avlKo(X)Q[($v  -Oy)«y 
+ «,  -8z)ez]e(X)dXd(£y  -Oy)d(t2-az)  (151) 
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By  setting  the  crystal  in  the  first  position  (parameters  /i01 , Q, , Q,  ) off  the  line  we 
obtain 


n(l,  l)c  = n(  1,  - l)c  = efXolw'1  ^R°0R°012d  cos J S$sd8zd8y 

+ R0R0l2d  cos^  | S\t8d^d?>y)  H57) 

and 

n(l»  *)C  = e(^o)w_1  iyRoRoi 2d  cosy?  Jstf&d8zd8y  +R%R%l2dcosV>  f sjfsd8zd8y) . 

(158) 

(Normally,  measurements  are  taken  on  each  side  of  the  line  and  averaged.)  Now  we  may 
obtain 


n(l,  DL  =n(l,  l)l+c  -n(l,  l)c 
n(l,  -1)L=  n(l,  -1  )^+c  - n(l,  -l)c 
n(l,i)L  =n(l,i)L+c  -n(l,i)c  , 
and,  from  Eqs.  (149)  and  (150), 

n,(l)L  = n,(l)L+c  - n,(l)c 
solely  due  to  line  radiation.*  We  now  let 


or 


uG  = 


con(l,  l)1,  con(l,  - 1 )L 


nt(l)L  n,(l)' 

?odo  I 


RlRoi  lStsQ?dS*dSy+RSR0i  KWM 6, 


U)G  — 


*01  sttQald8td6y*R*  j 


and  we  further  let 


H=n(  1,  i)L/n(l,  l)L  = n(l,  i)L/n(l,  -1)^ 


(159) 


(160) 


(161) 


(162) 


•A  note  of  caution:  The  measurements  here,  n(l,  1 )^+C,  n(l,  -l)L*C,  n(l,  i)L*C,  and  n,( l)L*C,  are 
functions  of  the  angle  between  the  centroid  of  the  incident  beam  and  the  crystal  planes,  great  care 
must  be  taken  to  ensure  that  this  angle  is  the  same  for  all  four  measurements. 
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or. 


H 


«i  J stS  ^id^d6y  + RoRoi  J s'l&Qnid^dK 

RoR°oi  |s#^d5A  + «i 


( 163) 


Now  source  functions  S^g  and  sjg'  will  have  the  same  angular  range,  which  can  be  made 
(and  usually  is)  much  larger  than  that  of  the  reflectivity  curve  of  a perfect  crystal,  the 
case  for  which  Q°  =£  Qn . We  therefore  let 


r/2  ol  L 
8 6 


Then 

g2  | Sl&  db;dby  = j sj'g  d8;d5y  (164) 

and  further, 

g2  Js^Q®d6*d5y-  | Q jd8id8y . (165) 

These  steps  would  have  been  impossible  if  we  bad  included  the  diffracted  continuum  as 
well  as  the  diffracted  line.  We  also  let 


/•2/?£  = fi£  (166) 

and  we  assume  that 

f2R°0X=R  Si  <167> 

i.e.,  the  two  crystals  are  the  same  in  this  respect.  (This  assumption  is  not  necessary  in 
cases  for  which  ^oi /^Oi  *s  known.)  We  may  now  write,  Kq.  (148), 

(2  + gl 

h = — 

1 + f2g2 

and 

n = f2n*g2) 

1 ♦ A V 

Elimination  of  produces 


(168) 


(169) 


/6  - 30/"*  + 3nf2  = 1 
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where 


1 +h  \ 
H )■ 


The  roots  of  this  equation  are 


and 


(171) 


Now,  if  f2  = 1 is  the  root,  H = h = 1.  Thus,  a = 1 and  we  find  that  all  three  roots  are 
the  same,  f2  = 1,  a trivial  result.  We  also  see  that  if  a < 1,  our  remaining  two  roots  are 
complex,  an  unphysical  result,  so  a > 1.  Comparison  with  our  diffraction  theories  shows 
that  the  only  reasonable  root  is 


r = 


3a  - 1 


3a-  lY 


]l/2 


-1 


(172) 


since  this  is  the  only  one  for  which  Rq/Rq  < 1,  in  general.  We  may  now  find  g2 , a 
measure  of  the  polarization  of  the  source.  Equation  (161)  may  now  be  written 


<joG  = R 


o 

0 


(1  +/V) 

(1  +f2g2) 


(1731 


from  which  Rq  may  be  found;  Rq  = f2R^  and  Rr0  = 1/2  ( R ^ + /?q)  readily  follow.  All 
the  integral  reflection  coefficients  are  thus  determined. 


5.2  Relative  Integral  Reflection  Coefficients 

By  following  the  procedure  in  Sec.  5.1  we  may  determine  the  integral  reflection 
coefficients  of  a crystal,  which  can  then  be  used  as  a standard.  When  such  a standard 
exists,  integral  reflection  coefficients  of  another  crystal  may  be  found  by  comparison. 

By  making  use  of  Kqs.  (143)  through  (146),  we  find  that  the  total  number  of  counts  in  a 
line  observed  by  our  detector  is  (for  C'  = f)) 


«(!)*'  = n{\),A(  -Atn,(l)(  =f(X0)w  1 RJ5(X0*  ‘S'/5 


+ «0<X0> 


(1741 
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and  (for  i//  = it/ 2), 

n(i)L  = n(i)L+c  - Atnt(i)c  = e(X0)«-*  Rn0(\ Q)  | Sjfrd8zd5y 

+ R°  Isjlfd6,d6y  . (175) 

Now,  by  adding  these  two  equations,  we  find 

2^n  = n(i)L  + n(l)L  = e(\0)u-'2RrQ(\0)  f sfgdSzdSy  (176) 

where  Sfg  = S^g  + S^g  . Thus  we  see  that  direct  comparison  between  the  measurements 
taken  with  the  standard  and  those  taken  with  our  sample  readily  provides  the  integral 
reflection  coefficient  Rq(Xq)  due  to  random  polarization  of  our  sample. 

Usually,  the  radiation  in  a spectral  line  is  not  greatly  polarized,  so  that  taking  the 
difference  of  n(i)L  and  nd)^  often  results  in  large  statistical  errors.  The  continuum 
radiation,  however,  may  be  polarized  to  a considerable  degree,  so  for  the  determination 
of  Rq  and  Rj  (or  rather,  of  their  difference)  we  will  use  the  continuum.  From  Eq.  (145) 
and  (146)  (if  A nct  = nt(i)c  - n,( l)c), 

Anf  =e(X0)2dcos^[R"(X0)-R£(X0)l  I (sj&  - sift  )dSzd8y . (177) 

Now,  direct  comparison  between  the  measurements  taken  with  the  standard  and  those 
taken  with  our  sample  provides  Rq(X0)  - Rq(X0).  By  using  Rq(X0)  as  determined,  we 
can  readily  find  Ry(\0)  and  Rq(X0). 

5.3  Reflectivity  Curve 

We  will  assume  here  that  we  have  an  unpolarized  source  of  line  radiation.  (A 
continuum  is  included  in  our  measurements;  but  we  have  seen  that  its  contribution  may 
readily  be  removed.)  In  front  of  this  source,  we  arrange  two  slits  of  equal  width  some 
distance  apart  and  uniformly  illuminated  by  the  source.  In  the  x-ray  region,  diffraction 
effects  due  to  macroscopic  slits  are  unimportant,  so  the  radiation  falling  on  a Bragg 
crystal  spectrometer  may  be  described  by  (assuming  random  polarization), 

sXi6  =*,(«,  )«(X  -X0)fc(Sy)  (178) 
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h(  6V)  = ° 


for  |6y  I > 50 


h(8y)  = j-  (80  +5y)  for  0 > 5y  > - 60 

h(6y)=  (S0  -5y)  forO<6y<60-  (17! 

We  will  require  (for  minimum  statistical  errors)  that  50  be  on  the  order  of,  or  less 
than,  the  diffraction  width.  From  the  detector  of  the  spectrometer,  we  obtain  a count 

rate, 


nt  = e(X0)  [ |V,(«,)cf«,]  |^(5y)K(X0)O°[0(Xo)  ' * - *>y 


■R£(X  0)Q7rl°(xo)-^-5yl}d5> 


(180) 


as  a function  of  crystal  angle  p.  Taking  the  second  derivative  of  n,(p)  with  respect  to  P 
and  integrating  by  parts  twice,  we  find  that 

— = const.  I'l8(6y  + 60)  - 25(6y)  + S(8y  - 6o)]Q[0(Xo)- P ~ «y  1 dSy  (181) 
dip2 

where  we  have  set  Q = R°0Q°  + The  integral  can  be  evaluated  trivially: 

d-^~  = const. {Q10(XO>  - P + Sol  -2QIO(X0)-  p]  Q[^(X0)  - V?  <182) 

dp2 

By  starting  with  data  far  from  the  line,  so  that  d2n,/dp 2 = o,  we  may  also  assume  that 

Q[0(Xo)  - p + 801  = Q10(XO)  - = Q|0(Xo)  - * ~5ol  = °- 

For  clarity,  we  set  7 = p - 0(Xq)  + xo-  Then  we  find 


Q(~7)  = 


d2nt(y-80) 


+ 2Q(-7+M-Q<-T'  + 26o>- 


Thus,  starting  at  a point  where  all  terms  equal  zero,  7 “ - 7o-  Q(>  fi0)  = 0,  we  find 

for  - 70  < 7 **  - To  + ^0> 


Q(-7)  = 


d2n,( 7 -?>0) 
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So  we  have  Q(70  ^ -7  > 7o  " 60)-  For  the  next  interval,  -70  + 60  < 7 < -70  + 2S0, 
we  find 

d2nt( 7 - 60) 

Q(-7)  = + 2Q(70  > - 7 + 60  > - 5o  )• 

d~rz 

For  the  next  and  all  succeeding  intervals,  we  must  use  all  the  terms,  since  none  are  cer- 
tain to  be  zero.  We  may  continue  on  in  the  same  way  until  we  have  the  complete 
shape  (except  for  a constant)  of  Q[0(X)  - ^] , the  integral  of  which  we  set  equal  to  unity. 
Because  of  statistical  considerations,  fitting  of  nt  may  be  required,  and  care  should  be 
taken  to  insure  that  the  second  derivatives  of  the  fitted  curve's  are  continuous. 


6.  INTEGRAL  REFLECTION  COEFFICIENTS: 

AN  EXAMPLE,  LITHIUM  FLUORIDE 

The  techniques  of  Sec.  5.1,  were  used  to  determine  the  first-order  integral  reflection 
coefficients  of  a lithium  fluoride  crystal  (200),  which  are  shown  here  as  an  example.  The 
detector  used  was  a flow  proportional  counter  filled  with  109£  methane  in  argon.  The 
pressure  to  which  the  counter  was  filled  became  the  largest  contribution  of  systematic 
errors.  Systematic  errors  were  determined  by  before-and-after  measurements  of  the  line 
intensity  diffracted  by  one  crystal.  The  errors  resulting  from  statistical  (Poisson)  varia- 
tions in  the  x-ray  source  were  taken  into  account  by  the  methods  described  in  Appendix 
A.  Errors,  both  systematic  and  statistical,  are  included  in  the  error  bars  of  the  figures. 

In  Fig.  8,  we  show  the  determined  ratio  Rq/Rq  as  a function  of  Bragg  angle  0 , 
together  with  that  expected  from  an  ideal  mosaic  crystal,  cos2 20,  and  that  expected  from 
an  ideal  perfect  crystal,  |cos20|.  Neither  of  these  idealizations  describes  this  ratio  over 
the  entire  wavelength  range  included  in  our  measurements.  In  Fig.  9,  the  three  integral 
reflection  coefficients  R g,  Rq,  and  Rr0  are  shown  as  functions  of  wavelength.  Since  all  the 
data  presented  in  Figs.  8 and  9 are  derived  from  the  same  set  of  measurements,  the  error 
in  each  datum  for  a particular  wavelength  is  not  completely  independent  of  the  error  in 
any  other.  For  example,  the  error  in  Rq(\  = 1.54  A)  is  not  found  by  taking  one-half  the 
square  root  of  the  sum  of  the  squares  of  /?q(X  = 1-54  A)  and  Rg(X  = 1.54  A). 

7.  CONCLUSIONS  AND  DISCUSSION 

In  the  preceding  sections,  we  have  discussed  the  usefulness  of  the  Bragg  single  crystal 
spectrometer  as  a tool  for  quantitatively  analyzing  the  spectra  of  x-ray  ('mission  sources. 
With  such  a device,  properly  calibrated,  we  may  determine  not  only  the  wavelengths  of 
spectral  features,  but  also  the  energy  contained  in  the  lines  and  continuum  as  well  as  the 
degree  of  polarization.  Line  shapes  can  also  be  determined  (with  some  difficulty), 
depending  on  the  intrinsic  line  width  and  spatial  extent  of  the  source.  The  Bragg  crystal 
spectrometer  is  therefore  a valuable  tool  for  diagnosing  x-ray  sources  (e.g.,  astrophysical 
and  laboratory  plasmas). 

We  have  presented  an  outline  of  the  theory  of  x-ray  diffraction  from  crystals  and 
discussed  the  qualitative  aspects.  Application  of  the  theory  to  real  crystals  (including 
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BRAGG  ANGLE.  6 


Fig.  8 — The  ratio  of  the  polarization  factor,  /<q/Hq,  fora  lithium  fluoride  crystal 
(200  cut).  Note  that  this  ratio  conforms  to  neither  the  perfect  crystal  diffraction 
theory  (|cos20|)  nor  the  ideally  imperfect  crystal  diffraction  theory  (cos22 0). 


Fig.  9 — First-order  integrated  reflection  coefficients  for  a lithium  fluoride  crystal  (200  cut ) as  a 
function  of  diffracted  wavelength.  Two  of  the  data  points  at  0.877  A are  displaced  in  wave 
length  for  clarity. 
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absorption  effects)  involves  extensive  numerical  calculation  even  if,  indeed,  the  theory 
applies  to  the  crystal  in  hand.  We  have  also  presented,  in  as  general  a form  as  possible, 
the  techniques  and  measurements  required  to  determine  the  characteristics  of  a single 
crystal  Bragg  spectrometer.  The  procedures  are  illustrated  using  a lithium  fluoride  crystal 
(200). 

j 
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Appendix  A 

STATISTICS  AND  PROPAGATION  OF  ERRORS 


Suppose  we  make  m independent  measurements.  Each  measurement  will  have  an 
associated  error  (e.g.,  if  N is  the  number  of  observed  photons  from  a source,  the  standard 
deviation  isy/N  if  the  photons  are  Poisson  distributed).  We  associate  orthogonal  unit 
vectors  e1(  e2,  e3  ...  em  with  the  m independent  measurements  Nx,  N2,  iV3,  ....  Nm. 
Denoting  the  standard  errors  by  olt  o2,  o 3,  ...,  am,  we  may  define  error  vectors  ax , o2, 
o3  ...,  om,  where  a,-  = o(e,-.  Now,  suppose  we  have  a calculation  based  directly  on 
measurements  of  AT,,  i.e., 


f = f(N1,N2,  ...  Nm ). 

Then,  approximately,  the  error  \<Jf\  in  f is  obtained  from 

_ V df 
°f=  L °‘- 

i=  1 ' 

Now,  suppose  we  calculate  an  additional  parameter: 

h-h{f.NltN2,...Nm  ). 
In  this  case  we  would  normally  take 


2=  f V.f02 

n t-*  \ df  m,  dNj  1 ' 
1=1  1 1 


(Al) 


(A2) 


(A3) 


( A 4 ) 


For  complicated  expressions  arrived  at  through  many  algebraic  steps,  the  propagation 
of  errors  in  such  a way  becomes  tedious  and  difficult.  The  same  result  is  obtained  if  we 
take 


oh  = 


'dh 


Of  + 


1*  l 


V’  dh 

dN:  °' 


(A5) 


and  let  ah  = \ah  I be  the  error  in  h.  Often,  the  final  results  desired  are  not  statistically 
independent  and  the  relationship  of  the  errors  in  the  final  result  might  also  be  required 
for  future  calculations.  Keeping  m-dimensional  error  vectors  is  often  not  practical.  By 
using  the  Schmidt  orthogonalization  procedure,  we  reduce  the  dimension  of  the  error 
vectors  to  at  most  the  number  of  final  results  desired. 
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